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ABSTRACT 

We report on the determination of the lithium abundance [A(Li)] of 52 solar-like stars. 
For 41 objects the A(Li) here presented corresponds to the first measurement. We have 
measured the equivalent widths of the 6708 A lithium feature in high-resolution spec¬ 
troscopic images (i? ~ 80 000), obtained at the Observatorio Astrofisico Guillermo 
Haro (Sonora, Mexico), as part of the first scientific observations of the revitalized 
Lunar and Planetary Laboratory (LPL) Echelle Spectrograph, now known as the 
Cananea High-resolution Spectrograph (CanHiS). Lithium abundances were derived 
with the Fortran code MOOG, using as fundamental input a set of atmospheric pa¬ 
rameters recently obtained by our group. 

With the help of an additional small sample with previous A(Li) determinations, we 
demonstrate that our lithium abundances are in agreement, to within uncertainties, 
with other works. Two target objects stand out from the rest of the sample. The star 
BD-l-47 3218 (Tes = 6050±52 K, A(Li) = 1.86± 0.07 dex) lies inside the so-called 
lithium desert in the the A(Li)-reff plane. The other object, BD-l-28 4515, has an 
A(Li) = 3.05±0.07 dex, which is the highest of our sample and compatible with the 
expected abundances of relatively young stars. 
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1 INTRODUCTION 


In convective stars, the photospheric lithium abundance 
A(LiLI decreases with time, because vertical motions trans¬ 
port Li atoms to deeper layers hot enough (« 2.5 x 10® K) 
to burn it by proton capture. 

The Sun has a very low pho tospheric A(Li) compare d 
to nearby solar analogs (e. g., iLambert fc Reddvl l2004h . 
and it is more than 2 dex lower than the meteoric abun- 
dance (A(Li)o = 1. 05 ± 0.10 vs. A(Li)n,et = 3.26 ± 0.05; 
lAsplnnd et al.l LOPOl ll. This difference cannot be explaine d 
by standard models (e. g., iD’Antona &i Mazzitellil Il994ll . 
consequently, several processes have been propo sed: for 


instance, mass loss dSwenson fc Fau lkner 1 992ll. angu¬ 


lar momentum loss and rot a tion Jschatzman fc Baelin 


1991 

: IPinsonneault et all 19921: Delivannis & Pinsonneault 

1997 

: lEggenberger. Maeder & Mevned 2010lL magnetic 

fields (Eggenberger. Maeder & Mevned 2010l: iLi et al. 

20141) or gravity waves JCarcia Looez & Soruit 

T^; 

Montalban & SchatzmarJl2000l: Charbonnel & Talon 

[2^, 


* E-mail: valdivia@inaoep.mx 

^ A(Li) = log(n(Li)/n(H)) -|- 12, where n is the number density 
of atoms. 


all of which induce an extra mixing. 

More recently, it has also been suggested that the pres- 
ence of planet s cou ld affect the evolut i on of Li abundance. 
iBouvieil (l2008ll and lEggenberger et all ll2012l ~l proposed that 
a long lasting star-disk interaction during the pre-main- 
sequence stage could slow down the rotation of the host star, 
enhanci ng the rotational mi x ing a nd thus destroying more 
lithium. iTheado fc Vauclaiil (120121 1 investigated the effects 
of disk matter accretion, that induce diffusive instabilities 
resulting in an extra mixing episode. ICastro et al.l ll2009l l 
pointed out that a transfer of angular momentum due to 
planetary migration could i ncrease the mi x ing eff iciency and 
Li depletion. Additionally, iTakeda et all ll2010ll confirmed 
that Li depletion in solar-like stars is more efficient when 
stellar rotation decelerates. 


Even though some authors fe. g.. Ilsraelian et all 

2004; 

Gonzalez! 20081: Israelian et all 20091: 

Gonzalea 

2014; 

Delgado Mena et all l2014l: iFigueira et al. 

I 2 OI 4 I: iGonzalez 

2015|1 claimed to have found that stars 

with exoplanets 


show an excess of Li depletion, the existence of a possible 
A(Li)- planet corre l ation is still matter o f debate, as different 
works llR vanll2()( H: iLuck fc Heite ill2006l: Ichezzi et^l2010l : 


iBaumann et all 2010l : Ramirez et all l2012h had not found 
such correlation. Increasing the number of main-sequence 
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(MS) stars with known Li abundance is therefore helpful to 
reach more robust statistical results. 

In this work, we present the determination of lithium 
abundances, measured from the 6708 A resonance line, for 
52 solar-like stars of intermediate brightness. For 41 of the 
targets our analysis represents the first Li measurement. The 
working sample and observations are described in Section 2. 
In Section 3, we detail the procedure followed for deriving 
abundances and, in Section 4, we discuss the results. 


2 SAMPLE DEFINITION AND 
OBSERVATIONS 


The stellar sample is somewhat inhomogeneous, but reflects 
a continuation of a global project aimed at determining stel¬ 
lar parameters (Teff / logg / [M/H]) as well as abundances 
of individual chemical species of a sample of Sun-like stars 
(types G0-G3) on the MS. As an ini tial step, we have se¬ 
lected 30 stars from our previous work (iLopez- Valdivia et al.l 
[2 o3 ): 12 objects of which correspond to metal-rich stars 
(SMR; [M/H]i^0.16 dex). This latter sub-sample is in fact 
composed by all the SMR stars that are sufficiently bright 
to allow high resolution observations with the instrument 
briefly described below. The sample is complemented with 
11 stars, that acc omplish the criteria in spect ral type and 
luminosity clas s of iLopez- Valdivia et al.l ((20l3), included in 
the catalog of ICasagrande et ah 1201 ill for which no pre¬ 
vious Li determination has been reported. For compari- 
son purposes, we als o observed 11 stars from the work of 
iRamirez et al.l ll2012l l who provide their Li abundances. 

The full sample is listed in Table [T] where we give 
some basic data and the results of our analysis. Golumns 
1-13 provide the star identification, the spectral type, the 


leadin g atmospheric parameters and their uncertaintie s 
(from Ihopez- Valdivia et al.l l2014l : iGasagrande et all I 2 OIIII , 
the 6708 A Li line equivalent width (EW) and its uncer¬ 
tainty, as described in the following section, and the last 
column gives a label for the source of stellar parameters. 

The spectroscopic data were collected at the 2.1 meter 
telescope of the Observatorio Astrofisico Guillermo Haro, lo¬ 
cated in Mexico, using the Cananea High-resolution Spectro¬ 
graph (CanHiSL as part of the first scientific observations of 
this instrumenl[j, at this facility. GanHiS is equipped with fil¬ 
ters, that provide access to ~40 A-wide wavelength intervals 
in single diffraction orders. Our observations were centered 
at 6710 A, with a spectral resolving power of i? ~ 80000 and 
a typical signal-to-noise ratio (S/N) of about 100. All the 
data collected (at least 3 images for each star) were reduced 
with IRAF, following standard procedures: bias substrac- 
tion, flat-field correction, cosmic-ray removal, wavelength 
calibration through an internal UNe lamp and, finally, con¬ 
tinuum normalization. We then transformed all the spec- 
tra to the rest fra me, using the spetrum of the Sun by 
iKurucz et al.l lll984| j as template. For each star we co-added 
its spectra weighted by their S/N. In Fig. [T] we present the 


^ CanHiS was donated in 2009 by the University of Arizona to 
our institute; it was originally known as LPL Echelle Spectro¬ 
graph and used at the 61-inch Catalina Telescope to carry out 
observations of Solar System bodies at spe cific spectral region s 
(for more details about its original design see iHunten et al ■I 1991 I) . 


spectra of the Li region of the 52 solar-like stars observed in 
this work. 


3 ABUNDANCES DETERMINATION 

In this Section we will present the procedure that we carried 
out for measuring the Li abundance of the stellar sample. 


3.1 Photospheric parameters and model 
atmospheres 


We collected the Tea, the surface gravity (log;?), the 
global metallicity ([M/H]) and their uncertainties (see Ta- 
ble m column 4 to 9 ) from two works. The first one is 
ILop^- V aldivia et al.l (I2014I L where Tea, logg and [M/H] 
were simultaneously determined using a set of Lick-like 
indices and intermediate-r esolution spectros c opic observa- 
tions. The second work is ICasagrande et al.l (1201 il l, where 
the Teff is obtained through the infrared flux method 
dBlackwell fc Shalh^l 19771 : ICasagrande et al1l2006ll , the log g 
from the fundamental relation, which involves bolometric lu¬ 
minosity and Teff and finally the metallicity as [Fe/H] from 
a calibration of Stromgren colour^. The last column of Ta¬ 
ble [ 1 ] reports the source of the stellar parameters for each 
star of our sample. 

By means of t he 22 obj ect in common i n the 
iLopez- Valdivia et al.l 1I2OI4I ') and ICasagrande et all (I2OIII1 
samples, we investigated possible systematic differences in 
the Tea scales, which is the most important parameter for Li 
abundance determination (see Sect. |T3T|. The data, plot¬ 


ted in Fig. (2] indeed show a tendency of ICasagrande et all 
ll201ll) temperat ures to be slightly h i gher, mostly at the 
cooler end. Since iLdnez- Valdivia et aH ll2014l ') Tea scale was 
constructed to optimi ze the parameters of G0- G3 MS stars 
and considering that ICasagrande et ahl ll201lll stated that 
their Tea are on average 100 K hotter that previ ous anal¬ 
yses, we trans formed the Casagrande et ^ (1201 ill temper¬ 
atures to the iLopez- Valdivia et all 1 2014h scale through a 
linear regression of data displayed in Fig. [2] 


t eff ,cor • 


Teff.Cas — 1247.86 
080 


( 1 ) 


where Teff.Cas is _ the temperature reported in 

ICasagrande et all (1201 il l and Teff,cor is the corrected 
one. 

We then deter mined a model atmospher e for each star 
by interpolating the lCastelli &: Kurucd (l2003l l grid of models 
in the Tefj, logg and [M/H] 3-dimensional space. For all 
cases, we used the microturbulence velocity ^ = 2 km s“^ 
(see also Section [3.3.1l) . 


3.2 Equivalent widths 

We measured the EW of the 6708 A lithium line with a 
flux fitting procedure. We simultaneously fit the Li and 
the adjacent Fe l line at 6707.44 A using three Gaussian 
functions: two of them are necessary to reproduce the two- 
component (simply called Lii and Li 2 ) asymmetric Li line, 

® We used the overa ll metallicity reported in the electronic ver¬ 
sion of the catalog of ICasagrande et al.l ll201lh . instead of [Fe/H]. 
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Table 1. Physical parameters and Li abundance of the stellar sample. Temperatures from ICasagrande et al.l (1201 il l were corrected following 
Eq.m 


Object 

V 

SType 

Teff 

(K) 

a 

(K) 

logff 

(dex) 

a 

(dex) 

[M/H] 

(dex) 

(7 

(dex) 

EW(Li) 

(mA) 

crEW(Li) 

(mA) 

A(Li) 

(dex) 

(TA(Li) 

(dex) 

S 

HD 5649 

8.70 

GOV 

5830 

52 

4.45 

0.22 

-0.08 

0.04 

27.5 

3.6 

2.05 

0.08 

L 

BD+60 402 

10.26 

GOV 

5985 

72 

4.30 

0.40 

0.22 

0.09 

58.0 

9.5 

2.58 

0.12 

L 

HD 13403 

7.00 

G3V 

5565 

204 

4.74 

- 

-0.29 

- 

< 5.5 

- 

< 1.04 

- 

C 

HD 16894 

8.02 

G2V 

5500 

70 

4.05 

0.30 

-0.10 

0.09 

6.4 

1.0 

1.1 

0.09 

L 

BD+60 600 

8.65 

GOV 

5655 

47 

3.95 

0.20 

0.20 

0.07 

78.0 

1.2 

2.51 

0.04 

L 

HD 232824 

9.52 

G2V 

5900 

67 

4.15 

0.35 

0.16 

0.08 

61.7 

5.7 

2.56 

0.08 

L 

HD 237200 

9.66 

GOV 

6045 

55 

4.25 

0.32 

0.18 

0.05 

18.3 

3.1 

2.03 

0.09 

L 

HD 26710 

7.18 

G2V 

5815 

47 

4.55 

0.20 

-0.04 

0.04 

60.8 

7.1 

2.47 

0.08 

L 

HD 31867 

8.05 

G2V 

5590 

57 

4.40 

0.25 

-0.10 

0.06 

< 4.9 

- 

< 1.03 

0.08 

L 

HD 33866 

7.87 

G2V 

5481 

123 

4.33 

- 

-0.07 

- 

11.0 

1.6 

1.28 

0.13 

C 

HD 41708 

8.03 

GOV 

5998 

58 

4.55 

- 

0.08 

- 

38.1 

1.2 

2.36 

0.05 

C 

HD 42807^ 

6.44 

G2V 

5617 

80 

4.53 

- 

-0.11 

- 

43.8 

5.5 

2.10 

0.10 

C 

HD 77730 

7.39 

G2V 

5698 

80 

4.13 

- 

-0.05 

- 

18.9 

2.8 

1.75 

0.10 

C 

HD 110882 

8.87 

GIV 

5880 

50 

4.40 

0.25 

-0.28 

0.04 

38.8 

0.8 

2.26 

0.05 

L 

HD 110884 

9.11 

G3V 

5905 

87 

4.30 

0.40 

-0.26 

0.08 

< 3.9 

- 

< 1.20 

- 

L 

HD 1115ld^ 

7.35 

GIV 

5723 

80 

4.31 

- 

0.12 

- 

11.4 

1.7 

1.53 

0.10 

C 

HD 111540 

9.54 

GIV 

5840 

47 

4.20 

0.25 

0.14 

0.05 

23.1 

5.2 

1.97 

0.12 

L 

HD 124019 

8.56 

G2V 

5685 

57 

4.65 

0.25 

-0.18 

0.06 

19.0 

2.6 

1.73 

0.08 

L 

HD 126991 

7.90 

G2V 

5360 

107 

3.15 

0.40 

-0.34 

0.14 

< 4.7 

- 

< 0.79 

- 

L 

HD 127913 

9.14 

G2V 

5475 

65 

4.45 

0.20 

-0.12 

0.08 

17.5 

1.2 

1.50 

0.07 

L 

HD 129357 

7.83 

G2V 

5775 

52 

4.30 

0.22 

-0.14 

0.05 

< 4.4 

- 

< 1.14 

- 

L 

HD 13094^ 

5.88 

F9IV-V 

5885 

80 

4.42 

- 

-0.09 

- 

93.8 

10.0 

2.84 

0.12 

C 

HD 13514^1 

8.35 

GOV 

5997 

80 

4.14 

- 

-0.02 

- 

20.3 

2.8 

2.04 

0.10 

C 

HD 135633 

8.46 

GOV 

6095 

67 

4.25 

0.40 

0.22 

0.06 

69.0 

0.9 

2.79 

0.06 

L 

HD 14038^1 

8.57 

G2V 

5735 

60 

4.60 

0.27 

-0.16 

0.08 

< 5.6 

- 

< 1.21 

- 

L 

HD 145404 

8.54 

GOV 

5920 

82 

4.43 

- 

-0.16 

- 

37.2 

4.6 

2.28 

0.10 

C 

HD 152264 

7.74 

GOV 

6177 

73 

4.09 

- 

0.02 

- 

53.6 

6.9 

2.70 

0.10 

C 

BD+29 2963 

8.42 

GOV 

5865 

55 

4.70 

0.22 

0.00 

0.04 

< 5.9 

- 

< 1.36 

- 

L 

HD 15696^1 

7.97 

GOV 

6105 

96 

4.42 

- 

-0.03 

- 

18.9 

1.4 

2.09 

0.09 

C 

HD 168874^ 

7.01 

G2IV 

5696 

80 

4.41 

- 

-0.05 

- 

65.4 

9.6 

2.41 

0.12 

C 

BD+28 3198 

8.66 

G2V 

5840 

35 

4.00 

0.17 

0.24 

0.05 

59.1 

8.1 

2.49 

0.09 

L 

HD 182407 

7.77 

GOV 

5953 

70 

3.94 

- 

0.25 

- 

< 3.6 

- 

< 1.22 

- 

C 

TYC 2655-3677-1 

9.93 

GOV 

6220 

47 

4.15 

0.27 

0.28 

0.05 

38.9 

6.2 

2.54 

0.10 

L 

HD 18723'P 

6.88 

G2IV-V 

5687 

62 

4.42 

- 

-0.04 

- 

89.2 

9.8 

2.63 

0.10 

C 

HD 187548 

7.98 

GOV 

6170 

42 

4.80 

0.17 

-0.24 

0.03 

49.3 

2.4 

2.63 

0.05 

L 

HD 333565 

8.75 

GIV 

5990 

52 

4.45 

0.27 

0.12 

0.05 

44.2 

6.5 

2.43 

0.09 

L 

HD 228356 

9.07 

GOV 

6055 

37 

4.00 

0.20 

0.16 

0.05 

63.4 

5.6 

2.71 

0.07 

L 

HD 193664^ 

5.93 

G3V 

5942 

112 

4.47 

- 

-0.11 

- 

30.2 

2.6 

2.19 

0.11 

C 

BD+47 3218 

8.70 

GOV 

6050 

52 

4.05 

0.30 

0.16 

0.06 

12.6 

1.30 

1.86 

0.07 

L 

HD 19996(P 

6.21 

GIV 

5895 

80 

4.24 

- 

0.22 

- 

61.0 

8.1 

2.54 

0.11 

C 

HD 201860 

8.65 

GOV 

5620 

107 

3.66 

- 

-0.02 

- 

30.9 

0.3 

1.93 

0.10 

C 

HD 210460 

6.19 

GOV 

5357 

80 

3.58 

- 

-0.17 

- 

< 4.8 

- 

< 0.80 

- 

c 

HD 210484^ 

7.59 

GIV 

5878 

95 

4.19 

- 

-0.01 

- 

19.8 

2.7 

1.93 

0.11 

c 

TYC 3973-1584-1 

10.74 

G2V 

6000 

45 

4.45 

0.22 

0.20 

0.05 

47.2 

6.4 

2.48 

0.09 

L 

TYC 3986-3381-1 

10.37 

G2V 

5855 

57 

4.15 

0.25 

0.26 

0.07 

79.1 

5.1 

2.68 

0.07 

L 

HD 239928 

8.69 

G2V 

5870 

67 

4.40 

0.30 

0.06 

0.06 

17.3 

2.4 

1.86 

0.09 

L 

HD 212809 

8.64 

G2V 

5975 

55 

4.55 

0.27 

0.16 

0.05 

33.5 

2.6 

2.27 

0.06 

L 

HD 217924 

7.22 

GOV 

5900 

52 

4.85 

0.20 

-0.10 

0.05 

17.2 

3.6 

1.87 

0.11 

L 

BD+28 4515 

8.73 

G2V 

5580 

40 

3.50 

0.17 

-0.22 

0.06 

146.0 

6.8 

3.05 

0.07 

L 

HD 218633 

8.19 

G2V 

5385 

70 

3.50 

0.30 

-0.32 

0.10 

< 6.4 

- 

< 0.96 

- 

L 

HD 218730 

7.32 

GOV 

5922 

83 

4.45 

- 

0.09 

- 

54.8 

9.5 

2.50 

0.12 

C 

HD 220334 

6.62 

G2V 

5801 

91 

4.13 

- 

0.14 

- 

87.9 

3.0 

2.72 

0.09 

C 


^ Star with previous determination of A(Li). 
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O 

+ 


n3 

i 2 

o 


1 1 1 1 1 1 1 1 1 1 1 1 
HD5649 ' ' 

1 1 1 1 1 1 1 1 1 1 1 1 
HD110882 ' ' 

1 1 1 1 1 1 1 1 1 1 1 1 
HD152264 ' ' 

1 1 1 1 1 1 1 1 1 1 1 1 
HD199960 ' ' 

BD+60402 1 1 

,y^ 

HD110884 1 1 

BD+292963 1 M 

HD201860 1 V . 


1 1 

1 1 


HD13403 1 If 

HD111513 1 1 

HD156968 | | 

HD210460 1 1 

1 1 

HD16894 

1 1 

HD111540 

1 1 

HD168874 

1 1 

HD210483 





BD+60600 1 1 

HD124019 1 1 

BD+283198 'W 

TYC3973-1584-1 ' ' , 

HD232824 I 1 / 

G166-31 1 1 

HD182407 1 H 

TYC3986-3381-1 I Y , 



- 


HD237200 . ^ i V 

1 1 

HD127913 , 1 

1 1 

TYC2655-3677-1 , , 

HD239928 | y 

1 1 

1 1 


1 1 

HD26710 

HD129357 | 

HD187237 ' ' 

HD212809 

^11 

HD31867 'V 

T 1 

HD130948 1 1 

HD187548 'W 

HD217924 1 ’ 

1 

HD33866 | | 

HD135145 |\|/ 

HD333565 | 1 

BD+284515 | | 

1 1 

11 



HD41708 I I 

HD135633 , . 

HD228356 j j 

HD218633 |||| 

HD42807 ' ' 

- 

HD140385 ' W 

HD193664 ' ' 

HD218730 ' V 

^ t. I I ... 


jy II 


HD77730 I I 

T 1 

HD145404 1 1 

BD+473218 1 1 

HD220334 1 V 

1 1 

1 1 

1 1 

Fe Li 

1 1 1 1 1 1 1 1 1 1 1 1 

1 1 

1 1 

1 1 

Fe Li 

1 1 1 1 1 1 1 1 1 1 1 1 

1 1 

1 1 

1 1 

1 1 

Fe Li 

1 1 1 1 1 1 1 1 1 1 1 1 

1 1 

Fe Li 

1 1 1 1 1 1 1 1 1 1 1 1 

6706 6708 

6706 6708 

6706 6708 

6706 6708 


Figure 1. Normalized spectra of our stellar 
is indicated by vertical dashed lines. 


wavelength (A) 

sample, centered around the Li line. The position of the Li line and of the adjacent Fe line 


while the third is needed to quantify the Fe contamination. 
The parameters of the Gaussian functions (central wave¬ 
length, width and amplitude) were modified within appro¬ 
priate values for each star to obtain the best fit, which was 
also validated through visual inspection. The integral of the 
sum of the Lii and Li 2 Gaussians provides the EW(Li). 

We adopted a Monte Carlo procedure to estimate the 
error in the measurement of EW(Li). For each star, we deter¬ 
mined the standard deviation of the flux in two wavelength 
regions, 6705.4-6707.0 A and 6708.6-6709.7 A, bracketing 
the Li line. We used this value as the width of a Gaussian 
distribution which provides the random noise value that we 
added to the observed spectrum and we repeated the fitting 
procedure previously described. We carried out 1000 realiza¬ 
tions to obtain the EW distributions of the three Gaussian 
functions, whose standard deviations we added in quadra¬ 


ture to obtain the error (TEW(Li). The EW(Li) values, along 
with their errors, are reported in columns 10 and 11 of Ta¬ 
ble [T] 


3.3 Lithium abundances computation 

To carry out the determination of the Li abundances, we 
made use of the driver abfin d of the Febru ary 2013 version 
of the Fortran code MOOG llSnederJll973l L which performs 
an adjustment of the abundances to match a single-line 
equivalent width. This code requires a model atmosphere 
and the physical parameters of a list of absorption lines to 
compute atomic abundances. We assumed the Li line to be 
a single electronic transition centered at A6707.84 A with 
log gf = -1-0.167, that corresponds to the average wave¬ 
length and the sum of the gf values of the two compo- 
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Fig ure 2. Comparison of T^ ff from [L6Dez-VaIdivia et alJ l|2014l') 
and ICasagrande et ^ ll201lf) for 22 solar-like stars in common. 
The dashed line shows the best linear regression, while the solid 
line indicates the one-to-one relation. 



A(Li) (This work) 

Figure 3. Co mparison of th e lithiu m abundance of the stars in 
common with iR'amfrez etU] (120121') . HD 187237 (see discussion 
in text) is represented by an empty square, while the solid line 
indicates the one-to-one relation. 


nents (Lii: A=6707.76 A, logff/=-0.009; Liz: A=6707.91 A, 
logg/=—0.309), obtained from the Vienna Atomic Line 
Database (VALD: lKupka et al.lll999l l. 


3.3.1 Error of Li abundances 

The main sources of error in the determination of the Li 
abundance are related with uncertainties on the stellar pa¬ 
rameters and on the EW measurement. In order to have 
an estimation of the errors introduced by the atmospheric 
parameters, we conducted some theoretical tests. We var¬ 
ied, one at a time, the four main stellar parameters (Teff, 
log^, [M/H] and ^), maintaining the other three fixed, and 
compute the Li abundance, assuming an EW(Li) = 50 mA, 
to measure the variations in A(Li). The chosen fixed values 
were 5750 K, 4.5 dex, 0.0 dex, and 2.0 km s“^, for Teff, logg, 
[M/H] and respectively. We obtained the following results: 
a change in effective temperature of ATes = 100 K produces 
a lithium abundance variation of AA(Li) = 0.09 dex, while 
a very small change in abundance of AA(Li) = 0.01 dex is 
obtained if we modify the other parameters by Alogp = 
0.5 dex, A[M/H] = 0.15 dex and A^ = 1 km s“^. These 
results, very similar to those from previous works (e. g., 
IChen et al.ll200l] : ICharbonnel fc PrimLaj[2005l i. indicate that 
the uncertanty in the Tea is by far the dominant source of 
error. Furthermore, our assumption of a single value for ^ 
of our models does not introduce a significant systematic 
error in the Li abundance. This choice is also supported , 
for instance, by the work of iGrav. Graham fc Hovu (l200lll . 
who derived a f between 1 and 2 km s“^ for their sample of 
F7-G2 MS stars. 

We calculated the contribution of the EW uncertanty 
to the error on the Li abundance [(TA(Li)] by computing the 
A(Li) for the values EW(Li)±(TEW(Li) and dividing their 
difference by half. The total CTA(Li) is the sum in quadrature 
of the errors on T^a and EW. All results are reported in 
Table [1] 


4 DISCUSSION 


In order to put the Li abundance of our sample in the con¬ 
text of the known results fo r solar-like stars, we used the 
work of iR.amfrez et ahl (|2012l i as a reference, as it is, to date, 
the largest homogeneous compilation of lithium abundance 
of MS and subgiants stars. They derived stellar parameters 
and Li abundances for 671 stars and included measurements 
from the literature to construct a catalog of 1381 FGK dwarf 
and subgiant stars. Analyzing the distribution of data in 
the A(Li)-refr plane, they confirmed that, in the interval 
5950 ^ Teff ^ 6100 K, stars appear to be neatly separated 
into two groups of high and low Li abundance, creating an 
avoidance regio n, called lith i um de sert, as was suggested for 
the first time bv IChen et al.l (1200 ill . These authors proposed 
that the low-Li side of the Af Lil-Teff diagram is populate d 
by evolved lithium dip objects (iBoessaard fc Tripiccdligsil l . 
that might bu r n the Li during the pre-MS stage. However, 
I Ramirez et al.l (l2012l l argues against this suggestion, based 
on their analyisis of t he mass-metallicity corr elation (see also 
iBalachandranI [l995l : iGummings et ah 2012h of their whole 
sample, which includes the IChen et al.| ll2001^ data. They 
also show that there are no young stars (age<2 Gyr) at the 
lower side of the lithium desert, suggesting that the mecha¬ 
nism that depletes surface lithium in the objects that popu¬ 
late that region s hould take place durin g the MS or subgiant 
phases. However. iRamfrez et al.l (l2012fl do not further spec¬ 
ulate on the nature of this mechanism. 

In total, we have 11 stars in common with 
IRamfrez et ah I (I2OI2I I: however, for the star HD 140385 we 
have only upper limits and was excluded from the compari¬ 
son. In Fig. |3] we present the comparison for the remaining 
common 10 stars, whose A(Li) are reported in Tabled The 
star HD 187237 presents the largest discrepancy of 0.45 dex. 
In order to achie v e the same lithium abundance reported 
in IRamfrez et al.l ll2012ll . this object should have half the 
EW measured in our work, which would represent a 4.bo- 
anomaly. Since, the EW measurement for this star is one of 
our best determinations, and the stellar parameters used in 
both works are not so different, we do not have a clear cause 
for the discrepancy. In the following analysis, we exclude this 
star too. In general, we found a good agreement between 
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Table 2. The lithium abundances and stellar parameters of the stars in common 


with lRami'rez et alJ <l2012l i. 


Object 

A(Li) 

Our 

(jA(Li) 

Our 

T'eff 

Ramirez 

log 9 

Ramirez 

[Fe/H] 

Ramirez 

A(Li) 

Ramirez 

(jA(Li) 

Ramirez 

HD 42807 

2.10 

0.10 

5674 

4.51 

-0.04 

1.92 

0.09 

HD 111513 

1.53 

0.10 

5811 

4.32 

0.11 

1.44 

0.03 

HD 130948 

2.84 

0.12 

5942 

4.39 

-0.11 

2.86 

0.07 

HD 135145 

2.04 

0.10 

5858 

4.09 

-0.08 

1.80 

0.04 

HD 156968 

2.09 

0.09 

5901 

4.36 

-0.14 

1.92 

0.02 

HD 168874 

2.41 

0.12 

5858 

4.44 

0.01 

2.68 

0.10 

HD 187237 

2.63 

0.10 

5792 

4.43 

0.02 

2.18 

0.10 

HD 193664 

2.19 

0.11 

5902 

4.42 

-0.14 

2.22 

0.02 

HD 199960 

2.54 

0.11 

5833 

4.20 

0.13 

2.43 

0.10 

HD 210483 

1.93 

0.11 

5834 

4.14 

-0.12 

1.97 

0.15 


the A(Li) dete r mined in this work and those determined in 
I Ramirez et alJ l|2ni2l ~): on average, we foun d a difference of 
+0.05 dex (our minus iRamirez et al. I l2012l data). This off¬ 
set is partly caused by adopting different stellar parameters 
(mainly different Tefi) and/or different atomic data. 

Figure |4] presents the A(Li) vs. Tes distribution 
of our results and their comparison with the complete 
I Ramirez et al.l (|2012l l compilation. It allows us to address 
three issues. 


Super metal-rich stars. 

Our stellar samp le inc ludes 12 SMR stars of 
iLopez- V aldivia et al.l ll 2014h f fM/H]>+0.16 dex) plus 

2 from Icasagrande et al.l (1201 il l, indicated with red sym¬ 
bols in Fig. |4] R is worthwhile mentioning that SMR 
stars have become a particularly interesting class, not only 
because they are important for modeling the chemical evo¬ 
lution of the Galaxy, but also for their pos itive correlation 
with the presence of giant planets (e. g.. lOon zaled fl9^: 
[Santos, Israelian fc MaTOijl2()0ll : iFischer fc Valentill2005ll 
The distribution of the Li abundances of our SMR stars 
is similar to that of the rest of our sample and it also re - 
sembles that of the SMR stars in iRamirez et al.l ll2012ll . 
All but three objects lie in the high abundance range of 
2.25 < A(Li) < 2.80 dex. In the case of the coolest object, 
BD+60 600, its A(Li)=2.51 dex is quite higher than the 
median value for its T^s class. On the warmer side, 3 stars 
show a lower than average Li abundance (A(Li)<2.10 dex): 
HD 237200, BD+47 3218, and HD 182407. The case of 
BD+47 3218 is discussed below in more detail. 


BD+47 3218: the star in the lithium desert. 


The star BD+47 3218 has Teff=6050±52 K and 

A(Li)=1.86±0.07 dex: these values locate it inside the 
lithium desert, at la in the case of Tefr and at 1.1a in the 
case of Li abundance. If we consider our poss ible small over- 
estima tion of the A(Li) with respect to the IRamirez et al.l 
l|2012ll scale, the star would move further inside the 
lithium desert. Encouraged by this finding, we additionally 
searched the literature for new A(Li) determinations for 


MS and subgiant s stars from 2012 to date. In total, we 
found 17 works (iMishenina et 2012 ; ICummings et all 


I 2 OI 2 I ; iLi et al.l |2012| ; Monaco et alJ 2012 ; IPace etlff 20121 : 


Schaeuble & King 20121: IXing et al. 20121: Xing & Xingl 

2 OI 2 I: Boesg 

aard et al. 20131: iFrancois et al.l 20131: 

Takeda et al.l 

2013; 

Gonzalez I20l4 Delgado Mena et al.l 

2014 lOnehag et al. 

I 2 OI 4 iDobrovolskas et al.l l2014l: 


Onehag et al.l 20 14 Gonzale3l2015l l , whose results provide 2 


more stars within the lithium desert boundaries : HD 44985 
(Teft=6004 K and A(Li) = 1.87 dex; lGonzale3 l2014ll and 
the star #584 40 of M4 (refi=5979 K and A(Li) = 1.82 dex; 
iMonaco et al.|[2012l L Keeping in mind that the A(Li) values 
come from different works, we can point out two possible 
conclusions regarding the existence of a lithium desert 


i) If there actually is a short-lived mechanism responsibl e 
for the lithium desert, as suggested bv iRamfrez et akl ll2012ll . 
the stars BD+47 3218, HD 44985 and #58440 in M4 repre¬ 
sent a good benchmark to investigate it. 
ii) The lithium desert is just a statistical product. 


The Li-rich star BD+28 4515. 

The star with the highest lithium abundance in our sam¬ 
ple is BD+28 4515, with A(Li)=3.05 dex. Comparing with 
the IRamirez et al.l 1 I 2 OI 2 II data, this value is the highest 
in a bin of temperature of 400 K, from 5350 to 5750 K. 
Wit hin this Teg ra n ge the re is a total of 228 objects, from 
the IRamirez et al.l (I 2 OI 2 I I and our samples, emphazising 
the peculiarity of BD+28 4515. From our literature revi¬ 
sion, we found 3 more stars in the same in terval with 
A(Li )>3.0 dex: HD 287927 (A(Li) =3.24 dex; Xing et al.l 
l2012ll . BD+07 582 IAILil=3. 26 dex:IXing et al.[l2012iL and 
HD 285840 (A(Li)=3.02 dex: IXing fc Xing||2012f) . They are 
all classified as weak line T Tauri stars (IXingl I 2 OI 0 I I. Us¬ 
ing the available data in the literature, we constructed the 
spectral energy distribution (SED) of BD+28 4515 and com- 
par ed it, in Fig. [5l wi t h its synthetic SED interpolated in 
the lCastelli fc Kurucd 1 I 2 OO 3 II grid, adopting the stellar pa¬ 
rameters reported in Table [T] Even though we notice the 
lack of the distinctive infrared excess of a classical T-Tauri 
star, the AKARI point, which might reflect an emission of 
the 10 +m silicate feature, appears to be in excess, how¬ 
ever at 2a significance only. We cannot rule out a possible 
weak line T-Tauri nature for BD+28 4515. Its high lithium 
abunda nce is compatible with t he values of the Pleiades (see 
Fig. 2 of iDelivannis et ^l2002lL which may be an indication 
of a relatively young age. 






















































































































Lithium abundance in a sample of solar-like stars 7 


2 - 


X 

CD 

g 1 

< 


0 - 


-1 


^ • • •• • — 
• • •• 


i ••• • • • '• 

• i’iV. li i I 

• Of«TP :• 1 iV ii, ^ 


• H ’ ^ i 

.,;. ‘ ‘ ‘. 

i*' ii 

i * * . I 


; ; 
i 


4500 


5000 5500 6000 6500 

effective temperature (K) 


7000 


Figure 4. Distribution of lithium abundance with Teff. Black filled circles indicate stars in our sample with [M/H] <0.16 dex, red squares 
show our SMR stars, and downward triangles mark 3(T upper limits. Grey dots and downward arrows represent iRamirez et al.l l|2012l 'l 
determinations and upper limit s. The polygon s hows the so called lithium desert. At the bottom right of the panel, we show the average 
errors for this work (black) and iRamfrez et~^ ] 1 I 2 OI 2 I) (gray). 
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Figure 5. Optical and infrared SED of BD+28 4515 . The sym¬ 
bols show the observed photometry, while the solid line represents 
the synthetic photospheric emission. Error bars are plotted, but 
they are in most cases much smaller than the symbol size. 
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